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We show that for a transiting exoplanet accompanied by a moon which also tran- 
sits, the absolute masses and radii of the star, planet and moon are determinable. 
For a planet-star system, it is well known that the density of the star is calculable 
from the lightcurve by manipulation of Kepler's Third Law. In an analogous way, the 
planetary density is calculable for a planet-moon system which transits a star, and 
thus the ratio-of-densities is known. By combining this ratio with the observed ratio- 
of-radii and the radial velocity measurements of the system, we show that the absolute 
dimensions of the star and planet are determinable. This means such systems could 
be used as calibrators of stellar evolution. The detection of dynamical effects, such as 
transit timing variations, allows the absolute mass of the moon to be determined as 
well, which may be combined with the radius to infer the satellite's composition. 
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1 INTRODUCTION 

Determining the absolute dimensions of a transiting planet 
system remains one of the most challenging problems 
confronting the exoplanetary scientist. The most widely 
adopted approach is to exploit Kepler's Third Law and 
the lig htcurve morphology to determ ine the stellar den- 
sity, p* l|Seager fc Mallen-Orn elas 2003). Combining p* with 
spectroscopic inferences of the surface gravity and metallic- 
ity allows us to produce stellar evolutionary tracks which 
may be used to est i mate the mass and radius of the star 
(e.g. lYi et all (|200ll ). lBaraffe et al.l i|l99Sl )). This approach is 
somewhat unsavoury as it is highly model dependent, where 
the model is that of stellar evolution. 

If the spectral lines from the transiting body are re- 
solvable, then the absolute dimensions may be determined 
as regularly done in the eclipsing binary community for 
double-lined, det a ched, eclipsing bin ary systems (DDEBs) 
l|Andersenl l|l99ll ); lTorres~ et al.l (|2009l ) and references there- 
in). This app roach was also rec ently used for the planet 
HD 209458b |Snellen et al.ll2010l ) but resulted in an error 
of greater than 20% on the stellar mass. Another method is 
to use interferometry to measure th e stellar radius bu t this 
is limited to nearby, bright targets (jBaines et al.|[2007l ). 

A more subtle oppor tunity for determi ning the absolute 
dimensions was offered bv lAgol et ai](|2005h . who considered 
dynamical interaction between two planets, where one tran- 
sited the host star. The presence of the other planet induces 



E-mail: d.kipping@ucl.ac.uk 



transit timing variations (TTV) onto the eclipsing body 
(see also lHolman fc Murray! l|2005l )). By measuring both the 
TTV and lightcurve derived stellar density, a unique solu- 
tion for the stellar mass is achievable. For a planet-moon 
system, perturbation of the planetary motion also gives rise 
to TTV (|Sartoretti fc Sdmeide"rlll999ri. as well as transit du- 
ration variations (TDV) (|Kippindl2009al ). 

In this letter, we show that the absolute masses and 
radii of the star and planet and the absolute radius of the 
moon are determinable for a transiting planet-moon system 
to an excellent approximation, even without measurable per- 
turbation to the planetary orbit. The detection of dynamical 
effects from the moon, such as TTV, allows for the absolute 
mass of the moon to be found too and a more accurate so- 
lution for the dimensions of all three bodies. It is expected 
that the Kepler Mission dBorucki et al. 20091) wil l be sensi- 
tive to moons of > 0.2MqT"( |Kipping et al. 2JD09c) and thus 
such systems could be detected in the near future. These 
systems therefore would obviate the need for stellar evolu- 
tionary models and potentially act as evolutionary calibra- 
tors. 



2 OBSERVABLES 

2.1 Star-Planet System 

We start by considering the simple case of a planet-star sys- 
tem. The lightcurve is completely described by just a few 
free parameters. These are the time of mid-transit, t m id, the 
impact parameter, b, the orbital period, P, the ratio-of-radii, 
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p — Rp/R* and the semi-major axis in units of the stellar 
radius, a/7?*. For eccentric systems we require the eccentric- 
ity, e, and the argument of periapsis, uj, as well. In total, we 
have 7 free parameters. We note that time of mid-transit is 
essentially just a phase term. 

Although b a nd a/R* share a strong correlation 
l|Carter et al.ll20oH ). typically a/R* is determinable to ~10% 
or better. For space based telescopes the precision can be 
a few percent, or better (|Kipping fc Bakosl l20ld ). If we 
know a/R* and P, we may exploit Kepler's Third Law 
to determine the stellar densi ty, as first pointed out by 
ISeager fc Mallen-Ornelasl (|2003T ). 



p* +p pp 



2,n{a/R, 
GP 2 



(1) 



A typical approach is to either assume p* 3> p 3 pp or 
iterate through the stellar evolution models a few times to 
improve the estimation of the p 3 pp term. 



2.2 Star-Planet-Moon System 

The lightcurve of a planet-moon system has similar free pa- 
rameters to that of the planet alone. We assume that in the 
frame of the barycentre of the planet-moon system, the or- 
bits of both bodies are unaffected by the star. Therefore, 
both the planet and moon orbit the barycentre in Keplerian 
orbits, which in turn orbits the star with Keplerian motion. 
This can be considered as a nested two-body simplification 
of the more general three-body problem. 

In this work, the only dynamical statement we make is 
that Kepler's Third Law is satisfied in these nested frames 
over the time-scale of the orbital period of the bodies. The 
validity of this approximation is discussed in the appendix, 
where we find that systems satisfying / 3 <C (3/2) are well- 
described by our approximation (where / is the semi-major 
axis of the moon around the planet in units of Hill radii). 
In practice, this means that all moons with / < 0.5 will be 
well-described here, w hich is known to be the bounding limit 
of prograde satellites IjDomingos et al.l l2006). 

In this set-up, the planetary parameters are unchanged, 
but should be given a subscript p to denote we are referring 
to the barycentre of the planet-moon system now. For exam- 
ple, P becomes Pp, i.e. the orbital period of the planet-moon 
barycentre around the host star: 

■ Pp: Orbital period of the planet-moon barycentre 
around the star 

■ ap/R*: Semi-major axis of the planet-moon barycentre 
around the star, in units of the stellar radius 

■ bp: Impact parameter of the planet-moon barycentre 
across the stellar face 

■ p: Ratio of the planetary radius to the stellar radius 

■ tmid.B: Time of mid-transit of the planet-moon barycen- 
tre across the stellar face 

■ ep: Eccentricity of the planet-moon barycentre around 
the star 

■ up: Argument of periapsis of the planet-moon barycen- 
tre around the star 

The satellite requires the same parameters except that 
it is understood we are referring to the orbit of the moon 
around the planet-moon barycentre at all times. The transit 



of a single body across the star is insensitive to the longi- 
tude of the ascending node but for two bodies the relative 
difference between the longitude of the ascending node be- 
comes relevant. We therefore include one extra angle, Sis to 
account for this differenc^E 

■ Ps: Orbital period of the moon around the planet- moon 
barycentre 

■ as / R* : Semi-major axis of the moon around the planet- 
moon barycentre in units of the stellar radius 

■ 6s: Impact parameter of the moon across the planetary 
face in the reference frame of the planet-moon barycentre. 

■ s: Ratio of the moon's radius to the stellar radius 

■ 4>s: Phase angle between the time of the planet- moon 
barycentre's transit across the stellar face and the moon's 
transit across the star. 

■ es: Eccentricity of the moon around the planet-moon 
barycentre 

■ uis: Argument of periapsis of the moon around the 
planet-moon barycentre 

■ Sis- Longitude of the ascending node of the moon 
around the planet-moon barycentre 

As before, we expect as/R* and bs to share a strong 
correlation. Nevertheless, it should be expected that a 
unique solution to both as / R* and Ps will be possible with 
sufficient quality photometry and quantity of transits. The 
investigation of parameter correlations and errors remains 
outside of the scope of this short letter. 

It is easy to see how as/R* would reveal itself in the 
lightcurve though. For a moon on a circular orbit, as/R* 
would be revealed by the temporal separation between the 
moon's mid-transit time and the planet's mid-transit time. 
In practice, one would fully model the transit to include 
both bodies and fit for as/R*- 

The other critical pa rameter, Ps, is also known to be 
measurable. For example, iKippinj (|2009al ) showed that the 
ratio of the TDV to TTV amplitude directly provides Ps- 
Another method is to use the time difference between ex- 
omoon transit features. For short period moons, this can 
be determined even with a single tr ansit in some cases, as 
pointed out bv lSato fc Asadal ^200$ ). 

2.3 Dynamical effects 

The planetary motion is perturbed by the presence of 
the moon giving rise to tran sit timing variations (TTV) 
jSartoretti fc Schneider[[ l999') and transit duration varia- 
tions (TDV) (k ipping 2009a, b). To completely describe the 
planetary motion, we still require on extra parameter, the 
mass ratio Msp — Ms/Mp, which is directly determinable 
from timing effects. With this ratio, and all of the previously 
mentioned parameters, the lightcurve of a planet-moon sys- 
tem may be completely described. 



3 DENSITIES 

Before we present our equations for d etermining the den - 
sities of all three bodies, we note that ISimon et alJ l|2007D 



The s subscript is used for the moon for consis tency with pre- 
vious works on exomoons e.g. ISimon et al" (2007) 
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discuss how the ratio of the exomoon-to-exoplanet density 
may be determined using a modified definition of transit 
time variation (TTV) which the authors labelled TTV P , or 
photometric TTV. However, the influence of the moon on 
the stellar and planetary derived densities was not consid- 
ered and the need for a customized definition of TTV is 
obviated in the equations presented in this work. Also, the 
authors did not discuss determining the absolute dimensions 
of the star or planet. 



3.1 Planetary density 

As we saw earlier, the ability to be measure cib/R* and P B 
for a transiting planet means that we can determine p* using 
Kepler's Third Law. In the same way, the ability to measure 
as/R* and Ps for the moon-planet system allows us to de- 
termine pp. No te that an explicit assumption throughout 
this work, as for lSeager fc Mallen-Ornelas! HoO.Si), is that all 
bodies are perfect spheres. 

One subtle difference is that the distance as is defined 
from the planet-moon barycentre to the star, but as is the 
moon to planet-moon barycentre separation. Kepler's Third 
Law is applicable to motion in the frame of reference of one 
of the two bodies, and not the barycentre. The planet-moon 
separation is therefore given by asp = as(l + Msp)- Given 
that as/R* and Msp are observables, then asp/R* is an 
observable too. Using Kepler's Third Law, we have: 



PS 



Ps 



M s 



(4/3)7^1 



PpMsp 



(5) 



Feeding this back in to the planetary density equation 
gives an equation for pp purely in terms of the observables. 



pp = 



37r(q s /E*) 3 (l + Msp) 



GPlp 



(6) 



It should be stressed at this stage that it is possible to 
provide an excellent estimation of pp even if no detectable 
TTV/TDV effects exist by setting Msp « 1- Indeed, the 
fact that no dynamical effects are detected actually requires 
Msp I- However, we continue using the unapproximated 
forms and present ps in terms of the observables only: 



Ps 



37r(a s /fl») 3 (l + Msp) 2 Msp 
GPSs 3 



And finally p* 



3tt 



(a B /R* 

p2 
r B 



(as/R,f(l + M S p) 
1 s 



(8) 



asp _ 1 (GM B Ps\ 1/3 
R* 



R* V 4tt 2 / 
3n(a S p/R*) 3 
GPIp 3 

For convenience, we replace asp with as- 



PP+{-) Ps = 



Pp+ { -) PS 



Mas/R*) 3 (l + M S p) 3 

gpIp 3 



(2) 



(3) 



The above expression is comparable to that used for 
deriving the stellar density in equation (1). However, (1) 
must be modified to account for the planet-moon system 
too: 



«B 

R, 



(G(M B +M t 



)P 2 b 



4tt 2 



1/3 



p* + p 3 pp + s 3 p s 



2,TT(a B /R«f 
GPl~ 



(4) 



Unfortunately, both equations (3) and (4) involve other 
densities and are not written purely in terms of the observ- 
ables. A resolution to this is to consider the satellite density 
and then return to these equations. 



4 ABSOLUTE DIMENSIONS 

Armed with p, p* and pp, it is possible to derive Mp* — 
M P /M*. 



Mp* = 



M P 



= p 3 ^ 
p* 



(9) 



In the case of Msp <S 1, we have shown earlier how pp 
is derivable even without detecting any TTV/TDV effects. 
Therefore, M P /M, is also derivable without TTV/TDV un- 
der the same approximation. However, continuing using the 
full unapproximated equation, it is possible to write the 
mass ratio of the planet-moon pair to the star as: 



Mb* = Mp» + Ms* 
P§(a B /R*) 3 



M 



P 2 B {a s /R,) 3 (1 + Msp) 



(10) 



At this point, we invoke information from a non- 
lightcurve origin. The radial velocity semi-amplitude con- 
tains information on the masses in the system and may be 
used to break the degeneracy one would otherwise be stuck 
with. Taking the standard equation for the semi-amplitude 
of the RV signal, and accounting for the fact the reflex mo- 
tion of the star is now due to the combined mass of the 
planet-moon pair, we have: 



3.2 Satellite density 

If we assume that one has measured Ms/Mp from the 
lightcurve, due to dynamical effects, it possible to evaluate 
the density of the satellite from the lightcurve alone. 



Kt = 



2-kGM'b sin 3 i s 



(M*+M S ) 2 P S (1- 4)3/2 



(11) 



Re-arranging for the stellar mass and replacing Mb 
M*Mb*, we obtain: 
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Table 1. Summary of derivables in two cases. Column 2 
shows the case where no dynamical effects are detectable (e.g. 
TTV/TDV) but the planet and moon transit the host star. Col- 
umn 3 shows the case where dynamical effects are detectable. 



Parameter 



No TTV/TDV 



TTV/TDV 



Af* Approximate solution 

R„ Approximate solution 

Mp Approximate solution 

Rp Approximate solution 
Mg Upper limit only 

Rs Approximate solution 



Accurate solution 
Accurate solution 
Accurate solution 
Accurate solution 
Accurate solution 
Accurate solution 



{1 + Mb*) J P B Kt(l-e B 



2 \3/2 



(12) 



i Bt sin is 

This may be re-expressed in terms of the observables 



only: 



(clb/R* 



(as/R*) 9 



p5 
B 



(1 



e 2 B f /2 Kl 



2nGs\n A i B 



{a B /R*fP 2 s - {a s /R*f{l + Msp?P 2 b 



(1 + M.sp) 9 
Similarly, the stellar radius may be written as: 



R* = 



(a B /R*) 2 VT 



(13) 



(14) 



2tt sin i B (as/R,) 3 (l + Msp) 3 PB 

Once the stellar mass and radius is known, the mass and 
radius of the planet are easily found. Once again, all of these 
quantities are derivable to an excellent approximation even 
without TTV/TDV by simply using A4sp "C 1 in equations 
(13) and (14). However, one can see that in the absence 
of a moon we have (a s /R*)/Ps -t 0/0 and (13) & (14) 
become meaningless, as expected. Given that i?« and s = 
Rs/R* have been determined, then the absolute radius of 
the moon is also known. In the absence of any TTV/TDV, 
we would be able to place at least an upper limit satellite 
mass, which would constrain the internal composition of the 
moon. However, if dynamical effects are detected then the 
more accurate solutions in equation (13) and (14) may be 
used and the absolute dimensions of all three bodies are 
determinable (Table 1). 



5 CONCLUSIONS 

We have shown that if the semi- major axis of a moon's orbit 
around an exoplanet and its corresponding orbital period are 
measurable, as is expected for a transiting planet-moon sys- 
tem, then it is possible to determine the densities of the star 
and planet. The detection of dynamical effects in the transit 
lightcurve (e.g. TTV/TDV effects) allows for an accurate 
determination of the various mass ratios and the density of 
the moon as well. However, an excellent approximate ex- 
pression for the mass ratio of the planet-moon pair to the 
star is still derivable even in the cases where Ms/Mp <C 1 
and no detectable TTV/TDV effects occur. 



Combining this mass ratio with the radial velocity data 
allows for an excellent approximate solution for the abso- 
lute masses and radii of the star and planet and the ab- 
solute radius of the satellite. The detection of dynamical 
effects allows for the more accurate solution for the abso- 
lute dimensions of all three bodies. The equations are not 
limited to strictly planet-moon cases but naturally are ap- 
plicable to binary-planets as well, should they be detected 
l|Podsiadlowski et al.|[2Qlol V 

Aside from such systems being interesting for harbour- 
ing a moon, they are also very valuable to validate the stellar 
evolution models, which are usually invoked for transiting 
planet systems. Planet-moon systems could therefore act as 
evolution calibrators, in the same way as double-lined, de- 
tached, eclipsing binaries (DDEBs) do. The solution also 
means that the absolute mass and radius of the planet and 
moon are resolvable. This is particularly useful for determin- 
ing the composition of the moon. An icy moon (which may 
have become a water- world after planetary migration) would 
have a significantly lower density than a rocky p lanet which 
may have been captured (|Grindrod et al.ll2007l ). Therefore, 
the density determination alone al lows for some a ssessment 
of the possible origin of th e moon dNamounill2010l l and pos- 
sibly even its habitability (|Williams et alj|l997h . 

The ability to measure the absolute mass and radius of 
the moon is extremely powerful. The moon is likely to be 
sub-Earth mass and modern radial velocity measurements 
would not be ab le to detect a planet of equivalent mass 
|Lovis et al.l I2OO8I ) . Therefore, moons measured in this way 
could become the objects of smallest measurable mass out- 
side of the Solar System. 
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2 3 (rs - r F ) 
-n s a s - 



2 3 r S 
Tlpdp- 



(15) 

\T S - rp\° ' ' |r S |° 
Where rs and rp are the position vectors of the satellite 
and planet relative to the star, tip is the mean-motion of the 
planet around the star and ap is the corresponding semi- 
major axis. Similarly, ns is the mean-motion of the moon 
around the planet and as is the corresponding semi-major 
axis. To translate to the reference frame of the planet, we 
may substitute: 



rs - rp 
-rp 



(16) 
(17) 



These vectors now describe the positions of the satellite 
and the star, respectively, relative to the planet. It follows 
that in this non-inertial frame where the planet is at rest, but 
the coordinate axes point in fixed directions, the equation 
of motion for the satellite becomes: 



2 3 r 2 3 

~ nsas M3 + n P a P 



(r' - r) 



(18) 



The term to the right of the addition sign represents 
the disturbance to Kepler's Third Law from a simple nested 
two-body approximation. The disturbing function can be 
written as: 



APPENDIX 

We will here explore the conditions under which our nested 
two-body formulation is a good approximation to the three- 
body problem, which can only be solved exactly through nu- 
merical techniques. The motion of the satellite, as the lowest- 
mass object in our hypothetical system, will be the body 
which experiences the greatest departure from the simple 
nested two-body approximation. In this letter, the deriva- 
tion has operated under the assumption that in the frame 
of reference of the planet, Kepler's Third Law applies for 
the motion of the moon, independent of the star. We are 
interested in the conditions under which the motion of the 
moon departs significantly from this statement (i.e. the star 
perturbs the satellite's orbit) and therefore our expressions 
are no longer valirfl 

For a planet-moon pair at infinite orbital separation, 
the validity our assumption can be understood intuitionally. 
As the planet-moon pair moves in, the disturbance due to 
the star becomes greater and our approximation will break 
down. We may follow the standar d treatment used in lunar 
theory (|Murrav fc Dermottl [1999) to calculate the distur- 
bance to our approximation. Let us select a non-rotating 
reference frame with the star at the rest at the origin, for 
which we may treat the frame as inertial since the mass 
of the star is much greater than the planet or moon. The 
satellite's equation of motion is given by: 



2 Orbital variations occurring on time-scales much greater than 
the orbital period, such as precession of the apsides or nodes, has 
a negligible impact on our derivation. 



a 2 3 
<±> = n P a P 
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Irpp 



|<J>| 

1*1 



2 

: ripCLp 



([1 + {as/a P )Y 



2n P as 



(19) 



Where have assumed a low eccentricity system and ex- 
panded to first order for as/ap <C 1. For the disturbing 
function to be small, we therefore require: 



2 3 r 

I r \ 3 



> 1*1 



n| 2np 



(20) 

Where the final line gives the condition under which our 
derivation is ultimately valid and is equivalent to: 



P| < Pp/2 



(21) 



It is possible to write Ps in terms of Pp and the moon's 
or bital distance in units of Hill radii, /, using the relation 
of lKippinel (1200931 ) that P s /Pp ~ \// 3 /3. Putting this to- 
gether yields: 



f < 3/2 



(22) 



Defining <C to indicate an order-of-magnitude differ- 
ence, i.e. a factor of 10, this constrains / ^ 0.531. We note 
that this distance is larger than that predicted as the maxi- 
mum stable separation for a prograde satellite of f — 0.4895 
(|Domingos et al.ll2006l ). However, a retrograde moon can be 
stable at up to / = 0.9309. Under such conditions, the equa- 
tions presented in this work would become invalid. 



